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LONG-TERM  GOALS 

The  long-tenn  goal  is  to  develop  a  nonhydrostatic,  parallel  ocean  simulation  tool  that  is  capable  of 
simulating  processes  on  a  wide  range  of  scales  through  use  of  accurate  numerical  methods  and  high- 
performance  computational  algorithms.  The  tool  will  be  applied  to  study  highly  nonlinear  internal 
waves  in  coastal  domains  to  develop  an  improved  understanding  of  mechanisms  that  govern  their 
generation,  propagation,  and  dissipation. 

OBJECTIVES 

The  primary  objective  is  to  enhance  the  capabilities  of  the  SUNTANS  model  through  development  of 
algorithms  to  study  multiscale  processes  in  estuaries  and  the  coastal  ocean.  This  involves  development 
of  1)  improved  momentum  and  scalar  advection  on  unstructured,  staggered  grids,  2)  accurate  and 
efficient  algorithms  for  solution  of  the  nonhydrostatic  pressure,  and  3)  adaptive  grid  capabiliites  with 
adaptive  mesh  refinement  and  model  nesting. 

APPROACH 

This  work  focuses  on  the  continued  development  of  SUNTANS  (Stanford  Unstructured 
Nonhydrostatic  Terrain- following  Adaptive  Navier-Stokes  Simulator),  a  free-surface,  nonhydrostatic, 
unstructured-grid,  parallel  coastal  ocean  and  estuary  simulation  tool  that  solves  the  Navier-Stokes 
equations  under  the  Boussinesq  approximation  (Fringer  et  ah,  2006).  The  fonnulation  is  based  on  the 
method  outlined  by  Casulli  and  Walters  (2000),  in  which  the  free-surface  and  vertical  diffusion  are 
discretized  with  the  theta  method  which  eliminates  the  Courant  condition  associated  with  fast  free- 
surface  waves  and  the  elevated  friction  tenn  associated  with  small  vertical  grid  spacing  at  the  free- 
surface  and  bottom  boundary.  For  flows  with  extensive  wetting  and  drying,  advection  of  momentum  is 
accomplished  with  the  semi-Lagrangian  advection  scheme  (Wang  et  al.  2011a),  which  ensures  stability 
in  the  presence  of  cells  that  fill  and  empty  with  the  tides.  Scalar  advection  is  accomplished  semi- 
implicitly  and  continuity  of  volume  and  mass  are  guaranteed  for  the  hydrostatic  solver,  following  Gross 
et  al.  (2002).  The  theta  method  for  the  free-surface  yields  a  two-dimensional  Poisson  equation,  and  the 
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nonhydrostatic  pressure  is  governed  by  a  three-dimensional  Poisson  equation.  These  are  both  solved 
with  the  preconditioned  conjugate  gradient  algorithm  with  Jacobi  and  block-Jacobi  preconditioning, 
respectively.  Because  the  nonhydrostatic  component  of  SUNTANS  is  essentially  a  correction  to  the 
hydrostatic  component,  SUNTANS  can  be  run  seamlessly  in  nonhydro  static  or  hydrostatic  modes. 
SUNTANS  is  written  in  the  C  programming  language,  and  the  message-passing  interface  (MPI)  is 
employed  for  use  in  a  distributed-memory  parallel  computing  environment.  SUNTANS  employs  the 
generalized  length  scale  approach  to  Reynolds-averaged  turbulence  modeling  (Wang  et  al.  2011b). 
The  SUNTANS  grid  employs  z-levels  in  the  vertical  and  is  unstructured  in  plan,  which  enables  the 
resolution  of  complex  coastlines  and  topographic  features.  Unstructured  grids  also  enable  the  use  of 
high  grid  resolution  in  regions  of  interest  while  coarsening  the  grid  in  regions  where  grid  resolution  is 
not  required,  thereby  significantly  reducing  computational  overhead. 

WORK  COMPLETED 

We  have  begun  the  development  of  a  nonhydrostatic  isopycnal-coordinate  model  and  have  developed 
methods  to  filter  unwanted  grid-scale  oscillations  associated  with  C-grid  discretizations  on 
unstructured  grids.  Stability  of  the  nonhydrostatic  algorithm  that  is  implemented  in  SUNTANS  has 
also  been  proven.  We  have  tested  SUNTANS  ability  to  predict  internal  wave  generation  at  an 
idealized  ridge  and  the  results  match  theoretical  predictions.  Below  we  focus  on  the  nonhydrostatic 
isopycnal  method  and  the  internal  wave  generation  results. 

RESULTS 

Nonhydrostatic  isopycnal-coordinate  model 

The  paper  by  Vitousek  and  Fringer  (2011)  indicates  that  accurate  resolution  of  internal  solitary  waves 
requires  a  horizontal  grid  spacing  that  is  less  than  the  mixed-layer  depth.  For  domains  like  the  South 
China  Sea,  this  requires  0(100  million)  grid  cells  just  to  begin  to  capture  the  longest  nonhydrostatic 
waves.  To  reduce  the  computational  expense  associated  with  accurately  simulating  nonhydrostatic 
solitary  waves,  we  have  begun  developing  a  nonhydrostatic  isopycnal  coordinate  model,  following  the 
success  and  efficiency  of  traditional  (hydrostatic)  isopycnal  coordinate  models  applied  to  the  South 
China  Sea  (Simmons  et  al.  2011). 

Like  any  nonhydrostatic  solver,  the  numerical  method  we  employ  requires  the  solution  of  an  elliptic 
equation  for  the  nonhydrostatic  pressure,  which  significantly  increases  the  computational  effort. 
However,  because  most  of  the  internal  wave  energy  is  contained  within  the  first  baroclinic  mode,  the 
isopycnal  approach  gains  efficiency  by  reducing  the  required  number  of  vertical  grid  points  from 
0(100)  grid  points  in  traditional  ocean  models  (as  in  Zhang  et  al.  201 1)  to  0(1),  thereby  reducing  the 
total  number  of  grid  cells  for  the  South  China  Sea  problem  to  a  more  tractable  value  of  0(1  million). 
Another  clear  advantage  of  moving  grid  or  isopycnal  models  is  that  they  reduce  or  eliminate  spurious 
diapycnal  mixing  (Griffies  et  al.  2000,  Koltakov  and  Fringer  2012).  Isopycnal  models  should  also 
eliminate  the  spurious  energy  loss  that  occurs  during  solitary  wave  formation  that  was  diagnosed  by 
Hodges  et  al.  (2006). 

Our  model  formulation  is  based  on  that  described  in  Casulli  (1997)  in  which  the  free  surface,  vertical 
diffusion,  and  isopycnal  layers  are  discretized  implicitly  with  the  theta  method.  To  accommodate  the 
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nonhydrostatic  pressure  in  our  model,  the  baroclinic  terms  are  discretized  explicitly  in  time  with  the 
third-order  Adams-Bashforth  method.  The  wetting  and  drying  of  isopycnal  layers  is  made  possible  by 
the  use  of  Multidimensional  Positive  Definite  Advection  Transport  Algorithms  (MPDATA; 
Smolarkiewicz  and  Margolin,  1998).  These  methods  ensure  positive  layer  heights  provided  the 
advective  Courant  number  is  less  than  unity.  The  model  is  currently  being  developed  for  structured 
grids,  however  the  approach  could  easily  be  extended  to  unstructured  grids  and  implemented  as  an 
option  in  SUNTANS. 

Adcroft  and  Hallberg  (2006)  outlined  an  approach  to  nonhydrostatic  modeling  in  isopycnal 
coordinates.  However,  they  ultimately  concluded  that  such  an  approach  was  “pointless  since  non¬ 
hydrostatic  effects  tend  to  be  associated  with  over  turning  (e.g.  Kelvin-Helmholtz  instabilities).” 
However,  we  believe  that  our  approach  is  essential  for  accurate  resolution  of  nonlinear  and 
nonhydrostatic  internal  solitary  waves.  Although  the  model  is  in  its  early  stages  of  development,  it  is 
already  capable  of  forming  well-resolved  and  physically  realistic  nonlinear  internal  solitary  waves  in 
idealized  test  cases,  as  shown  in  Figure  1. 
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Figure  1:  Comparison  of  the  evolution  of  a  solitary-like  wave  using  the  hydrostatic  isopycnal  model 
(red)  and  the  nonhydrostatic  isopycnal  model  (blue)  with  the  horizontal  grid  resolution  equal  to  the 
upper  mixed-layer  depth  (similar  to  the  test  case  from  Vitousek  and  Fringer  2011).  The 
nonhydrostatic  isopycnal  model  forms  physically  realistic  solitary  waves  while  the  hydrostatic  model 

forms  waves  that  are  too  large  and  narrow. 
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Nonhydrostatic  and  nonlinear  internal  wave  generation 

We  have  developed  a  code  based  on  the  formulation  of  Kang  and  Fringer  (2012)  to  calculate  internal 
wave  energy  radiation  from  a  model  ridge  in  the  SUNTANS  model.  The  results  show  the  variation  of 
the  linear  internal  wave  energy  flux  as  a  function  of  the  tidal  excursion  parameter  (s2=tidal 
excursion/topographic  scale)  for  three  different  slope  regimes  based  on  the  value  of  the  criticality 
parameter  (si  =Topographic  slope/internal  wave  beam  slope).  Because  linear  theory  is  restricted  to 
small  tidal  excursions  (linear  flow),  small  topographic  slopes  (subcritical  slope),  finite  water  depth,  and 
simple  topography  shapes  (e.g.  Witch  of  Agnesi  profile;  St.  Laurent  et  al.  2003),  we  compare  the 
SUNTANS  results  to  those  from  the  Green’s  function  approach  which  extends  linear  theory  to 
topographies  of  arbitrary  shapes  and  slopes  (Petrelis,  2005;  Echeverri  and  Peacock,  2011).  The  iTides 
model  (Echeverri  and  Peacock,  2011)  is  used  to  compute  theoretical  predictions  of  internal  wave 
energy  flux  for  comparison  to  the  SUNTANS  model. 

To  compare  the  results  of  SUNTANS  to  those  of  the  iTides  model  (Echeverri  &  Peacock,  2011),  all 
tenns  in  the  SUNTANS  model  are  linearized  except  for  the  bottom  boundary  condition.  A  semidiurnal 
baro tropic  tide  is  forced  over  a  Gaussian  ridge  with  a  height  of  ho= 50  m  and  width  of  12  km  in  a 
constant  stratification  such  that  the  slope  of  the  internal  tidal  beams  is  1.6°.  The  SUNTANS  domain  is 
two-dimensional  with  a  length  L  =  1400  km  and  depth  H  =  1000  m,  with  a  horizontal  grid  resolution  of 
2  km,  a  vertical  grid  resolution  of  8.3  m  (120  vertical  grid  cells),  and  a  time  step  size  of  120  s.  As 
shown  in  Figures  2-4,  excellent  agreement  is  obtained  between  SUNTANS  and  iTides  for  various 
values  of  the  slope  and  excursion  parameters.  The  parameters  shown  in  Figures  3  and  4  are  the  same  as 
those  in  Figure  2  except  that  in  the  critical  case  H=  3500  m  and  ho  =  205  m  and  in  the  supercritical 
case  H  =  4500  m  and  ho  =  300  m.  The  SUNTANS  model  results  match  those  of  iTides  to  within  1%. 

As  shown  in  Figure  5,  when  the  nonlinear  and  nonhydrostatic  terms  are  included  in  the  SUNTANS 
model,  SUNTANS  produces  energy  fluxes  that  are  smaller  than  those  of  iTides.  We  are  analyzing  the 
terms  in  the  energy  flux  budget  to  determine  the  physical  reasoning  behind  this  discrepancy. 

IMPACT/APPLICATIONS 

High-resolution  simulations  using  nonhydro  static  models  like  SUNTANS  are  crucial  for  understanding 
multiscale  processes  that  are  unresolved,  and  hence  parameterized,  in  larger-scale  ocean  models. 
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Figure  2:  Comparison  of  internal  tide  energy  flux  using  hydrostatic  and  linear  SUNTANS  (solid 
line)  and  Hides  (dashed  line)  for  a  subcritical  Gaussian  ridge. 
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Figure  3:  Comparison  of  internal  tide  energy  flux  using  hydrostatic  and  linear  SUNTANS  (solid 
line)  and  iTides  (dashed  line)  for  a  critical  Gaussian  ridge. 
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Figure  4:  Comparison  of  internal  tide  energy  flux  using  hydrostatic  and  linear  SUNTANS  (solid 
line)  and  Hides  (dashed  line)  for  a  supercritical  Gaussian  ridge. 
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Figure  5:  Comparison  of  internal  tide  energy  flux  using  nonhydrostatic  and  nonlinear  SUNTANS 
(solid  line)  and  iTides  (dashed  line)  for  a  critical  Gaussian  ridge. 


6 


REFERENCES 


Adcroft,  A.,  and  R.  Hallberg,  2006,  On  methods  for  solving  the  oceanic  equations  of  motion  in 
generalized  vertical  coordinates,  Ocean  Modelling,  1 1(1),  224-233. 

Casulli,  V.,  1997,  Numerical  simulation  of  three  dimensional  free  surface  flow  in  isopycnal 
coordinates.  International  J.  for  Numerical  Methods  in  Fluids,  25(6),  645-658. 

Echeverri,  P.,  Peacock,  T.,  2010  Internal  tide  generation  by  arbitrary  two-dimensional  topography,  J. 
Fluid  Mech.,  659,  247-266. 

Griffies,  S.  M.,  Boning,  C.,  Bryan,  F.  O.,  Chassignet,  E.  P.,  Gerdes,  R.,  Hasumi,  H.,  and  D.  Webb, 
2000,  Developments  in  ocean  climate  modelling.  Ocean  Modelling,  2(3),  123-192. 

Hodges,  B.  R.,  Laval,  B.,  B.  M.  Wadzuk,  2006,  Numerical  error  assessment  and  a  temporal  horizon  for 
internal  waves  in  a  hydrostatic  model,  Ocean  Modelling,  13(1),  44-64. 

Kang,  D.  and  Fringer,  O.  B.  2012  Energetics  of  barotropic  and  baroclinic  tides  in  the  Monterey  Bay 
area,  J.  Phys.  Oceanogr.,  42  (2),  272-290 

Koltakov,  S.  and  O.  B.  Fringer,  2012,  Moving  grid  method  for  numerical  simulation  of  stratified  flows, 
International  J.  for  Numerical  Methods  in  Fluids,  published  online  29  August  2012, 
doi:  10.2002/fld.3724. 

Petrelis,  F.,  Llewellyn  Smith,  S.,  Young,  W.  R.  2006  Tidal  conversion  at  a  submarine  ridge,  J.  Phys. 
Oceanogr., 36  (6),  1053-1071. 

Simmons,  H.,  Chang,  M.  H.,  Chang,  Y.  T.,  Chao,  S.  Y.,  Fringer,  O.,  Jackson,  C.  R.,  and  D.  S.  Ko, 
2011,  Modeling  and  prediction  of  internal  waves  in  the  South  China  Sea.  Oceanography- 
Oceanography  Society,  24(4),  88. 

Smolarkiewicz,  P.  K.,  and  L.  G.  Margolin,  1998,  MPDATA:  A  finite-difference  solver  for  geophysical 
flows,  J.  Computational  Physics,  140(2),  459-480. 

Vitousek,  S.,  and  O.  B.  Fringer,  2011,  Physical  vs.  numerical  dispersion  in  nonhydrostatic  ocean 
modeling,  Ocean  Modelling,  40  (1),  72-86. 

Zhang,  Z.,  Fringer,  O.  B.,  and  S.  R.  Ramp,  2011,  Three-dimensional,  nonhydrostatic  numerical 
simulation  of  nonlinear  internal  wave  generation  and  propagation  in  the  South  China  Sea,  J. 
Geophys.  Res.,  116,  C05022. 

PUBLICATIONS 

Carter,  G.S.,  Fringer,  O.B.,  and  E.  D.  Zaron,  2012,  Regional  models  of  internal  tides,  Oceanography, 

25  (2),  56-65,  doi:10.5670/oceanog.2012.42. 


7 


Giddings,  S.N.,  Fong,  D.A.,  Monismith,  S.G.,  Chickadel,  C.C.,  Edwards,  K.A.,  Plant,  W.J.,  Wang,  B., 
Fringer,  O.B.,  Horner-Devine,  A.R.,  and  A.T.  Jessup,  2012,  Frontogenesis  and  frontal  progression 
of  a  trapping-generated  estuarine  convergence  front  and  its  influence  on  mixing  and  stratification, 
Estuaries  and  Coasts,  35  (2),  665-681,  doi:  10. 1007/sl2237-01 1-9453-z. 

Kang,  D.,  and  O.  B.  Fringer,  2012,  Energetics  of  barotropic  and  baroclinic  tides  in  the  Monterey  Bay 
area,  J.  Phys.  Oceanogr.,  42  (2),  272-290. 

Koltakov,  S.,  and  O.  B.  Fringer,  2012,  Moving  grid  method  for  numerical  simulation  of  stratified 
flows,  Int.  J.  Numer.  Methods  Fluids,  published  online  29  August  2012,  doi:  10.2002/fld. 3724. 

Koltakov,  S.,  Iaccarrino,  G.,  and  O.  B.  Fringer,  2012,  Inferring  bottom  bathymetry  from  free-surface 
flow  features,  SIAM  Conference  on  Uncertainty  Quantification,  Raleigh,  North  Carolina. 

Simmons,  H.,  M.-H.  Chang,  Y.-T.  Chang,  S.-Y.  Chao,  O.  Fringer,  C.R.  Jackson,  and  D.S.  Ko.  2011. 
Modeling  and  prediction  of  internal  waves  in  the  South  China  Sea.  Oceanography  24(4):88-99, 
doi:10.5670/oceanog.20i  1.97. 

Trigo  Cabrita  Gil,  G.,  and  O.  B.  Fringer,  2012,  Lagrangian-  and  Eulerian-mean  effects  in  progressive 
internal  gravity  waves,  2012  Ocean  Sciences  Meeting,  Salt  Lake  City,  Abstract  10029. 

Venayagamoorthy,  S.K.,  and  O.  B.  Fringer,  2012,  Examining  breaking  internal  waves  on  a  shelf  slope 
using  numerical  simulations,  Oceanography  25  (2),  132-139,  doi:10.5670/oceanog.2012.48. 

Vitousek,  S.,  and  O.  B.  Fringer,  2012,  Stability  and  consistency  of  nonhydrostatic  free-surface  models 
using  the  semi-implicit  theta-method,  Int.  J.  Numer.  Methods  Fluids,  in  press. 

Vitousek,  S.,  and  O.  B.  Fringer,  2012,  Grid  resolution  requirements  in  modeling  internal  waves,  2012 
Ocean  Sciences  Meeting,  Salt  Lake  City,  Abstract  12494. 

Walter,  R.K.,  Woodson,  C.  B.,  Arthur,  R.  S.,  Fringer,  O.B.,  and  S.  G.  Monismith,  2012,  Nearshore 
internal  bores  and  turbulent  mixing  in  southern  Monterey  Bay",  J.  Geophys.  Res.,  117,  C07017, 
doi:10.1029/2012JC0081 15. 

Wolfram,  P.,  and  O.  Fringer,  2012,  Mitigating  divergence-error  oscillations  in  triangular  C  grids  for 
nonlinear  and  nonhydro  static  flows,  Workshop  on  Multiscale  Modelling  of  Coastal,  Shelf  and 
Ocean  Dynamics  (IMUM),  Delft,  Netherlands. 


8 


